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ABSTRACT 

We study the structure and substructure of halos obtained in N— body simulations for a A Cold 
Dark Matter (ACDM) cosmology with non-Gaussian initial conditions (NGICs). The initial statistics 
are lognormal in the gravitational potential field with positive (LN P ) and negative (LN„) skewness; the 
sign of the skewness is conserved by the density field, and the power spectrum is the same for all the 
• simulations. Our aim is not to test a given non-Gaussian statistics, but to explore the generic effect of 

| positive- and negative-skew statistics on halo properties. From our low-resolution simulations, we find 

£NJ . that LN p (LN„) halos are systematically more (less) concentrated than their Gaussian counterparts. 

^sO ' This result is confirmed by our Milky Way- and cluster-sized halos resimulated with high-resolution. In 

addition, they show inner density profiles that depend on the statistics: the innermost slopes of LN p 
(LN„) halos are steeper (shallower) than those obtained from the corresponding Gaussian halos. A 
subhalo population embedded in LN p halos is more susceptible to destruction than its counterpart inside 
Gaussian halos. On the other hand, subhalos in LN„ halos tend to survive longer than subhalos in 
Gaussian halos. The spin parameter probability distribution of LN p (LN„) halos is skewed to smaller 
(larger) values with respect to the Gaussian case. Our results show how the statistics of the primordial 
density field can influence some halo properties, opening this the possibility to constrain, although 
indirectly, the primordial statistics at small scale. 



1. INTRODUCTION paper, we are interested in particular in the effects of NG- 

The inflationary-motivated Cold Dark Matter (CDM) ^8 on the density profiles of galaxy- and cluster-sized 

, , ./ ii i u j c u CDM halos as well as on the substructure abundance m 

cosmology describes rather well a large body of observa- ^ c i Vd [ 0S 

■ tional data at large scales and has served as an useful the- „ „ . . r ,. , 

;* , 4--if l-c - 4--i- i 4? j following, we present a short review of theoretical 

oretical framework for investigating galaxy formation and , .. , ' . , . , , , , 

, . . / r , • tv. i onno I?- • studies and observational constraints related to non- 
evolution (see tor recent reviews e.g., r/renk 2002; r lrmam „ . . . ,. . „ 
c a -i T3 r, nno \ t j-u i~itm\? i j-i, ■ i Gaussiamty, as well as a brief discussion of previous nu- 
& A vila- Reese 2003). In the CDM cosmology, the initial • , , , , , r . ■ , , . — , 
, ., a . ,. i . , ,. . ° r i i . merical works on large-scale structure formation using N G- 
density fluctuation field is the starting point for calculat- ^ ° ° 

ing the structure formation in the universe. The simplest 

assumption for the statistical distribution of the primor- 1.1. Generation of Non- Gaussian fluctuations 
dial density fluctuations is the Gaussian one, which is valid 

only in the limit of zero fluctuation amplitude; for finite Inthe contact of inflationary models, non-Gaussian fluc- 

rms fluctuations, the Gaussian statistics assigns a non-zero tu f 10ns w ° uld , be the f S natur ^ of non-lmeanties during 

probability to regions of negative density. In the context Ration. For the simplest single scalar field models the 

of simple inflationary models, the Gaussian distribution slow-roll conditions prevent non-Gaussian primordial fluc- 

may arise in a natural fashion. However, as we will see Rations from developing significantly (Kofman & Lmdc 

below, deviations from Gaussianity are possible, in partic- 19 , 87 ; Ortolan, Lucchin & Matarrese 1988, Matarrese, Or- 

ular for more complicated models of inflation, while from tolan & ^f™^ ^ arrow & Co es 1990; but see Ac " 

the observational point of view the results regarding non- ^ uavi T a e * aL 2002 ' who suggest that non-Gaussiamty 

Gaussianity are controversial. Moreover, it may happen m ^ d f elo P ft second order calculation m cosmological 

that the statistics depend on scale, so that measurements perturbations) . These constraints can be partially relaxed 

made at large scales are not necessarily representative of when special non-lmear conditions such as a non-vacuum 

the small scale statistics. An interesting question to ex- J nl ] 1™°^ t° U - ( ^ ^ azu f° H akel " 

plore is how non-Gaussian initial conditions (NGICs) do lar f dou 2000 > Gan S u 1 1 ' Martm & Sakellariadou 2002) or 

affect cosmic structure formation at small scales. In this a . deP arture from mvanance m the mflaton poten- 

tial (Lesgourgues, Polarski & Starobmsky 1998; Salopek 
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1992) are introduced. Deviations from a scale invariant 
spectrum indeed can make the primordial non-Gaussianity 
non-negligible (Acquaviva et al.). There are some pieces 
of evidence that the measured matter power spectrum 
presents a peculiar behavior: a sharp peak or break at 
k w O.OSh^Mpc (e.g., Einasto et al. 1997; Peacock 1997; 
Gaztahaga & Baugh 1998), an abrupt bending of the spec- 
tral index n s at scales smaller than k w 0.1h -1 Mpc (Pea- 
cock 1997), and/or a runing n s from n s > 1 at large scales 
to n s < 1 at small scales (Peiris et al. 2003). 

However, significant deviations from Gaussianity are 
more probable to occur in multi-field inflationary models. 
In this case, the deviations are produced by fluctuations 
of the field that does not drive inflation, and it is conse- 
quently not constrained by a flat potential. Generically, 
one gets a mixture of adiabatic and isocurvature perturba- 
tions, quadratic in a Gaussian field (e.g., Allen, Grinstein, 
& Wise 1987; Peebles 1999; Antoniadis, Mazur, & Mot- 
tola 1997; Linde & Mukhanov 1997). In some proposals, 
the mixing between adiabatic and isocurvature modes may 
effectively transfer the non-Gaussian features from one sec- 
tor to the other, rendering them more significant (Gordon 
et al. 2001, and references therein; Bernardeau & Uzan 
2002; Bartolo, Matarrese & Riotto 2002). An overview 
of the characteristics and limitations of the non-Gaussian 
contributions in the different inflationary models, either 
with a single or with several scalar fields is presented by 
Bernardeau & Uzan. 

A recently proposed alternative scenario consists in adi- 
abatic density perturbations which originate after inflation 
from the perturbation of a field that is not the inflaton, and 
does not dominate the dynamics, the so called curvaton 
(Lyth & Wands 2002; Lyth, Ungarelli & Wands 2003; see 
also Linde & Mukhanov 1997; Moroi & Takahashi 2001, 
2002). The perturbation is initially of the isocurvature 
type but, if the curvaton is long-lived, its energy density 
will have a relative increase with respect to radiation, and 
will finally generate a curvature perturbation. The curva- 
ton scenario can generate significant non-Gaussianity since 
the curvature perturbation is proportional to the pertur- 
bation in the curvaton energy density, which is in turn a 
linear combination of the curvaton-field perturbation and 
its square. In this model, large isocurvature perturbations, 
correlated with adiabatic perturbations, may arise in dark 
matter and baryonic matter. 

In contrast, in the topological-dcfect model, strong non- 
Gaussian (skew-positive) fluctuations in the density field 
rather than in the potential field may develop (e.g., Turok 
& Spergel 1990; Avelino, Shellard, & Wu 1998). Although 
this model has proved to be inconsistent with observations 
(e.g., Pen, Seljak & Turok 1997; Durrer, Kunz & Mel- 
chiorri 1998; Allen et al. 1997), a mixed model of inflation 
plus certain defects (cosmic strings and any other kind of 
global defects) seems to be still allowed by observations 
(Bouchet et al. 2002; Sakellariadou 2002 and references 
therein). 

1.2. Observational constraints 

As it has been just discussed, the theoretical prediction 
of the exact amplitude and precise form of non-Gaussianity 
in inflation is still an open question. Thus, we should turn 
to observations to look for its signatures and constraints 
(see Kamionkowski 2003 for a recent review on the sub- 



ject). Care must be taken in drawing conclusions about 
the Gaussian nature of the density field from observations; 
the uncertainty in the measurements and the limited cov- 
erage of present observations can hide cosmological non- 
Gaussian features. We should also be aware that the grav- 
itational evolution and some systematic effects introduce 
non-Gaussian signatures in an intrinsically Gaussian dis- 
tribution. We also would like to stress that there are no 
a priori strong arguments to postulate that the statistics 
of the primordial fluctuations should be the same at all 
scales. 

The most direct way to search for primordial non- 
Gaussianity (at large scales) is the Cosmic Microwave 
Background Radiation (CMBR). Many tests applied to 
the CMBR data resulted nearly consistent with Gaussian- 
ity (see e.g., Kogut et al. 1996, and Bromley & Tcgmark 
2000, and references therein, for COBE DMR maps; Po- 
lenta et al. 2002 for BOOMERanG maps; Santos et al. 
2002, and Wu et al. 2001 for MAXIMA maps; Park et 
al. 2001 for QMAP and Saskatoon maps; Komatsu, et al. 
2003 for WMAP results). However, non-Gaussian statis- 
tics have no generic signatures on the sky and different 
tests and data analysis, at different scales, may be better or 
worse suited for different types of non-Gaussian statistics. 
Therefore, it is important to consider carefully any new re- 
sult that supports the possibility of non-Gaussian features 
using different techniques and indicators (e.g., Ferreira, 
Magueijo & Gorski 1998; Magueijo 2000; Pando, Valls- 
Gabaud & Fang 1998; Chiang et al. 2003). Notice that 
Banday, Zaroubi, & Gorski (2000) have shown that the 
non-Gaussian signal detected in the COBE data by Fer- 
reira et al. was probably due to a systematic period in 
the COBE satellite flight — the "eclipse period" . See also 
Barreiro et al. (2000) for some comments on the other 
mentioned works. 

The difficulty of predicting the precise form of non- 
Gaussianity from inflation or other seed mechanisms, led 
to define and work with some parametric indicators of 
the deviation from Gaussianity in a generic model instead 
of trying to look for a particular statistics. An often 
employed indicator is the non-linear coupling parameter 
/nl that parameterizes the leading order non-linear/non- 
Gaussian correction to the primordial gravitational poten- 
tial: <& = 4>l + /nl[0|— < <fit >]' wnere 0l is a Gaus- 
sian random field. A similar definition is introduced for 
non-Gaussianity in the density field (see e.g., Verde et 
al. 2000). The recent first-year WMAP observations 
(processed with two statistics: the angular bispectrum and 
Minkowski functionals) have constrained /nl in the range 
-58 < / NL < 134 (95%) (Komatsu, et al. 2003). As a ref- 
erence, let us recall that in some non-standard inflationary 
models, as the multi-field ones, /nl can be as big as ~ 20 
(Verde 2001), while in the single-field models with some 
non-linearities, |/nl| <5 1 (Acquaviva et al. 2002). Com- 
bining WMAP data with ACBAR, 2dF, and HST results, 
strong constraints can be obtained in specific cases of cor- 
related primordial and adiabatic density fluctuations, such 
as curvaton models (Gordon & Lewis 2002). 

Tegmark, de Oliveira-Costa & Hamilton (2003) have 
performed an independent foreground analysis from the 
WMAP data. Their foreground-cleaned and Wiener- 
filtered maps were used by Chiang et al. (2003) for a 
non-Gaussianity test. These authors implement a phase- 
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mapping technique (Chiang, Coles & Naselsky 2002), 
which has the advantage of testing non-Gaussianity at 
separate multipole bands (angular scales). A multipole 
band which is considerable non-Gaussian could have an 
insignificant non-Gaussian contribution to the whole map 
and produce thus an overall Gaussian realization. The 
authors show that the foreground-cleaned map is against 
Gaussianity at all 4 multipole bands centered around 
£ = 150,290,400, and 500. The Wiener-filtered map, on 
the other hand, is Gaussian for I < 250 (k, 180h _1 Mpc) 
but non-Gaussian for i > 350 (<J 50h _1 Mpc). As it was 
said above, the statistics might change with the scale. 

The effects of primordial non-Gaussianity can be also 
studied in the distribution of mass in the Universe (in- 
termediate scales). The gravitational evolution not only 
amplifies the primordial density fluctuation field but also 
makes it non-Gaussian due to clustering. However, the 
non-Gaussianity generated by gravitational evolution can 
be calculated and subtracted from the present statistics 
(e.g., Verde et al. 2000). Some of the tests of non- 
Gaussianity include the abundance and distribution of 
galaxy clusters and galaxies, the evolution of these abun- 
dances, and the size-temperature relation of galaxy clus- 
ters. 

Since galaxy clusters form at the highest density peaks, 
which are most sensitive to deviations from Gaussianity, 
their abundance, and in particular the evolution of their 
abundance, can be used to probe the statistics. A den- 
sity distribution with positive (negative) skewness leads 
to more (less) high density peaks. In this field, a useful 
indicator of the amount of deviation from Gaussianity is 
the probability of obtaining a peak of height 3a or higher 
(at a given "filtered" scale) for a certain statistics relative 
to the Gaussian one: 



G, 



2tt 



f 3 °° p (y)dy 



(1) 



Robinson, Gawiscr & Silk (2000), and Koyama, Soda 
& Taruya (1999) (see also Chiu, Ostrikcr & Staruss 
1998; Willick 2000) constrain primordial non-Gaussianity 
through the evolution of cluster abundance. They use a 
Press-Schechter formalism extended to NGICs for predict- 
ing cluster abundances at different redshifts and for dif- 
ferent cosmological models. Koyama et al., using data 
sets at z = and z ~ 0.6 for cluster abundance and the 
correlation length from the APM survey, obtained the up- 
per bound fi m < 0.5 and the lower bounds G3 > 2 and 
og > 0.7 for a x 2 non-Gaussian distribution; in particu- 
lar for S7 m = 0.3 and £!a = 0.7, non-Gaussianity of the 
order of G3 ~ 4 (skew positive) is favored. Robinson et 
al. (2000), from several samples of clusters at redshifts 
between 0.05 and 0.8, find that for TO = 0.3, Gaussian- 
ity is consistent with the data, but a wide range of non- 
Gaussian lognormal models also fit the data, with values 
of G3 < 4 (G3 < 6 when including cosmological constant) 
acceptable at the 2tr level. It must be stressed that current 
constraints are not conclusive due to the difficulty in the 
determination of the mass of high-redshift clusters. 

A long non-Gaussian tail would give not only an en- 
hanced cluster abundance, but also a larger scatter in the 
formation redshifts and sizes of clusters with a given mass. 
A skew-positive distribution implies that halos collapse 



during a broader range of redshifts than in the Gaussian 
case, leading to a larger scatter of clusters properties, such 
as size and temperature. Verde et al. (2001a) find that the 
predicted scatter for Gaussian initial conditions is consis- 
tent with the observed one and constrain G3 to be <; 4. 

Galaxies at high redshift are another kind of rare objects 
in the universe, so their abundance is considered for testing 
non-Gaussianity by Verde et al. (2001b). Results are still 
very preliminar and uncertain. On the other hand, the bis- 
pectrum of the galaxy distribution from 2dF (Verde et al. 
2002) and from SDSS (Szapudi et al. 2002) surveys have 
been found consistent with Gaussian initial conditions. 

Finally, the statistics of large-scale galaxy clustering can 
also be used to constrain non-Gaussianity. Most of the 
previous cosmological numerical simulations with NGICs 
were aimed at studying precisely galaxy clustering. In the 
next subsection we review some of these works. In the 
present paper, we will introduce a new potential test for 
cosmological models with NGICs at the smallest scales, by 
exploring the inner structure and substructure of the dark 
matter halos. 

1.3. Previous numerical works 

In an attempt to overcome the problem of lack of power 
at large-scales of the standard CDM cosmology, among 
other alternatives, models with NGICs were proposed and 
tested by means of N-body simulations. Moscardini et al. 
(1991, hereafter MMLM; see also Messina et al 1990) con- 
structed certain multiplicative non-Gaussian statistics for 
the peculiar gravitational potential, and used them as ini- 
tial conditions of A— body cosmological simulations. The 
outcome of the simulations were analyzed by different sta- 
tistical tests at large-scales and compared with observa- 
tions. It was found that both the clustering dynamics 
and the present day texture are mostly sensitive to the 
sign of the initial skewness of density fluctuations. Pos- 
itive models cluster to a lumpy structure with small co- 
herence length, while negative models build up a cellular 
structure with large coherence length, large voids, and en- 
hanced peculiar motions. In a subsequent paper, the same 
authors showed that some relevant differences may appear 
for models with the same sign of the skewness but dif- 
ferent statistics (Coles et al. 1993). For example, they 
showed that their CDM skew-negative \ 2 model agrees 
rather well with the 2D topology inferred from the Lick 
Catalogue, while the skew-negative lognormal model (as 
well as the skew-positive models) give discrepant results. 
From the analysis at smaller scales, Matarrese et al. (1991) 
found that groups in skew-negative models are preferen- 
tially loose systems, slightly biased over the background, 
with high velocity dispersion and clear signal alignment 
with neighboring groups. Groups in skew-positive mod- 
els are highly biased compact systems with lower velocity 
dispersion. In general, in the series of papers mentioned 
above, the authors concluded that standard CDM skew- 
negative models (predominance of primordial underdense 
regions) retain the advantages of the Gaussian biased stan- 
dard CDM counterpart, and improve on several of their 
problems. 

An extensive numerical study of the influence of NG- 
ICs on large-scale structure formation has been done also 
by Weinberg & Cole (1992). These authors attempted to 
isolate the effects of NGICs from other factors by con- 
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sidering a range of power-law power spectra and different 
cosmological models. Similar to Matarrese et al. (1991), 
they found that skew-positive models form structure by ac- 
creting mass onto rare high peaks, leaving large volumes 
almost unperturbed, while skew-negative models develop 
structures rich in expanding voids, sheets and filaments. 
The authors conclude that skew-positive models produce 
unrealistic large-scale features, while skew-negative mod- 
els create attractive structure but, for the low-density uni- 
verse (O m = 0.2) cluster velocity dispersions are too low 
and the voids are probably excesively large. 

ft should be noted that since the time these works ap- 
peared (early 90's), deep progresses in the determination 
of the cosmological parameters, the observational data on 
galaxy clustering, and the computational capability were 
made. Time is ripe to explore the effects of NGICs on 
stucture formation by means of high-resolution N— body 
simulations for the preferred (concordance) cosmological 
model. This paper presents by first time a numerical study 
about the structure of halos formed under NGICs. 

In §2 we present the non-Gaussian models to be used 
here, while in §3 the implementation of the NGICs in the 
N— body code is described and the different low- and high- 
resolution runs are presented. In §4.1, the halo concentra- 
tions and density profiles for the non-Gaussian and Gaus- 
sian runs are presented and compared among them. The 
spin parameter distribution of the low-resolution halos and 
the angular momentum distribution of the high-resolution 
halos is analyzed in §4.2. The ellipticity of the latter halos 
is analyzed in §4.3. In §4.4, results on the subhalo popu- 
lation in the galaxy and cluster halos are presented. The 
summary and the discussion are given in §5. 

2. NON-GAUSSIAN MODELS 

Since there is no preference for any particular non- 
Gaussian statistics (§1.1), we chose arbitrary non- 
Gaussian models. Our main interest is to study generically 
the influence of skew-positive and -negative NGICs on halo 
properties. To get the initial conditions for our simula- 
tions, we have followed closely MMLM, who constructed 
the non-Gaussianity in the gravitational potential fluctu- 
ation field, <&, instead of the density field. The potential 
fluctuation $ is related to the linear mass fluctuation via 
the Poisson equation. We first generate a Gaussian re- 
alization in Fourier space of a w random field with the 
primordial (scale-invariant) power spectrum 



P w (k) = Bk~ 



(2) 



where k is the comoving wave number, and B is the nor- 
malization factor. The zero-mean Gaussian field w is ob- 
tained in configuration space by inverse- Fourier transform- 
ing w(k). The non-Gaussian field for the perturbations in 
the potential, <&(a;), is then built as the convolution of a 
real function t(x) with a stationary, zero-mean random 
field ci(x), which is obtained from a non-linear transfor- 
mation of the w(x) Gaussian field. For a lognormal statis- 
tics the relation between w and the primordial potential 
fluctuation <ft is given by 



4>{x) = A 



MX) 



(3) 



transfer function T(k) used to get the processed power 
spectrum at the recombination epoch. The Fourier trans- 
formed function $(£;), which is simply given by the prod- 
uct of <j)(k) with T 2 (fc), is subsequently processed by the 
Zel'dovich approximation algorithm. 

Different lognormal statistics are obtained depending on 
the amplitude of the w Gaussian field. Although MMLM 
set (w 2 ) = 1, we prefer to assume a variance slightly 
smaller than one because it reproduces better the theoreti- 
cal ACDM-dcnsity power spectrum. There is still freedom 
in choosing the sign of the constant A in eq. (3), with 
the sign determining the direction of the skewness. Our 
Gaussian realizations are also implemented using a </> field 
but in this case 



(x) = Aw(x). 



(4) 



We notice that because Gaussian models with oposite signs 
give rise to different realizations in a statistical sense, 
we had to simulate for each lognormal skew-positive LN p 
or lognormal skew-negative LN„ model the corresponding 
gaussian G p or G„ model. We also notice that because in 
general a statistical distribution for the gravitational po- 
tential does not imply the same distribution for the density 
field (MMLM), we should not expect the statistics of the 
density field to be lognormal, nevertheless the sign of the 
skewness is preserved. 

In order to compute the probability distribution func- 
tion (PDF) for the different statistics at the initial time 
Zi = 50 (after the application of the Zel'dovich approxi- 
mation), our 100h _1 Mpc box simulation (see below) was 
divided in 64 3 cells. The density in each cell was then 
computed using the cloud-in-ccll scheme. The skewness 



parameter £ 



3/2 



is 0.52, 3.90, and -1.27 for 



where A is a normalization factor, related to B in eq. 
(2). The Fourier transformed of the function t(x) is the 



the Gaussian, LN p and LN„ statistics, respectively, while 
the G3 parameter defined in eq. (1) is 3.05 for LN p and 
0.067 for LN„. Note that at z t the PDF of the density 
field of a Gaussian realization is not Gaussian; in fact, we 
found that the PDF could be fitted by a lognormal, skew- 
positive, distribution: non-linear modes appear as a result 
of the application of the Zel'dovich approximation. This 
is why the skewness of the density field in the Gaussian 
scenario is not zero anymore. 

3. NUMERICAL SIMULATIONS 

The Adaptive Refinement Tree (ART) TV-body code 
(Kravtsov, Klypin, & Khokhlov 1997) in its multiple-mass 
version was used to run a series of simulations aimed 
mainly at studying the structure of halos under different 
statistics. The ART code achieves high spatial resolution 
by refining the base uniform grid in all high-density regions 
with an automated refinement algorithm. The multiple- 
mass scheme is used to increase the mass and spatial res- 
olution in few selected halos and its implementation is 
described in detail elsewhere (Klypin et al. 2001). All 
A— body simulations were performed in a low-density flat 
ACDM model with the following parameters: Q m — 0.3, 
= 0.7, and h = 0.7. The power spectrum was normal- 
ized to erg = 1.0, where a% is the rms mass fluctuation in 
spheres of radius 8h _1 Mpc (the top-hat window is used). 

As a test of our method for generating NGICs, we first 
ran two simulations under G p and LN p statistics in a box 
of 20 h _1 Mpc for each side and 64 3 particles. Halos from 
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Fig. 1. — Projected particle positions in slices of one tenth the size of the 100 h _1 Mpc box at z = 0. Each slice contain about 
26,000 dark matter particles. We note that the particle distributions are compared at the same cosmological time and not when the 
clustering properties of the different scenarios are comparable, as has often been the case in the literature (e.g., MMLM). Although, 
the differences in the large-scale pattern between a lognormal realization and its Gaussian counterpart are not as striking as if 
they were compared at the same clustering epoch, we do see some differences: voids are emptier and more defined in the LN n than 
in their Gaussian scenario. On the contrary, in the LN P scenario voids are not as developed (they contain more DM particles) as 
in the Gaussian one. 



these simulations were also used for further analysis. We 
also ran two simulations with low-mass resolution in boxes 
of 12.5 and 100 h _1 Mpc on a side. We selected, for resim- 
ulation with high-mass resolution and for each statistics, 
two galaxy-sized halos from the 12.5h _1 Mpc box and one 
cluster-sized halo from the 100h _1 Mpc box. In an attempt 
to study galaxy-sized halos with different mass assembly 
histories (MAHs), we simulated with high resolution two 
halos, Gl and G2, with different concentrations since it 
is known that the concentration of the halo is related to 
its MAH (Avila-Reese et al. 1998, 1999; Wechsler et al. 
2002). We notice, however, that in a box as small as the 
one we chose for the search of galaxy-sized halos, only 
two or three halos can be found to satisfy simultaneously 
the following criteria: (a) comparable mass, (b) extreme 
concentrations, (c) isolation, and (d) with a counterpart 



in the corresponding lognormal realization. Halos G1q p 
and G2q p , for example, obey roughly these rules. They 
have masses that differ by a factor of three and NFW ( 
Navarro, Frenk & White 1997; Eke, Navarro & Steinmetz 
2001) concentrations, c v ; r , of 16.4 and 11.4, for G1g p and 
G2q p , respectively. The reason why we did not find a 
galaxy-sized halo with a very low concentration is because 
very often these halos are found accompanied with halos of 
comparable mass. In the rare case in which we do find an 
isolated halo with a low c v ; r , it turns out to be too small 
or to lack a "partner" in the lognormal statistics. 

The bound density maxima (BDM) group finding algo- 
rithm (Klypin & Holtzman 1997) was used to locate halos 
and subhalos in all simulations. The BDM algorithm iden- 
tifies positions of local maxima in the density field at the 
scale of interest and applies physically motivated criteria 
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Table 1 
Simulation Parameters 



statistics 


Lbox 


m p 


Resolution 




Vmax 


c l/5 


A' 


Halo name tag 




(h^Mpc) 


(h" 1 M ) 


(h^kpc) 


(h^M©) 


(kms* 1 ) 


(io- 2 ) 




G p 


100.0 


5.0 x 10 9 


3.0 


6.2 x 10 14 


1407 


6.1 


2.82 


C1 GP 




100.0 


5.0 x 10 9 


1.5 


8.3 x 10 14 


1643 


7.2 


1.37 


C1ln p 




i nn n 


c; n v 1 n9 

O.U X 1U 


O.U 


1.0 X 1U 




D.o 


*3 70 
o. t J 




LN n 


100.0 


5.0 x 10 9 


6.1 


1.3 x 10 14 


781 


5.2 


1.22 


ClLNn 


Gp 


12.5 
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to find out whether a group of particles forms a gravita- 
tionally bound halo. 

An overview of some of the characteristics of the stud- 
ied halos are shown in Table 1. The symbol with which we 
denote the statistics is shown in column (1). The mass of 
the least massive particle is shown in column (3). It is only 
a function of the size of the box (col. [2]) because the cos- 
mological model and the number of mass species are fixed. 
Column (4) shows the formal force resolution, measured 
by the size of a cell in the finest mesh. In general, the 
lognormal skew-positive (skew-negative) simulations have 
a higher (lower) resolution than the Gaussian counterparts 
because simulated halos in the former case reach a higher 
(lower) central density than in the latter one. The mass 
M„ within the virial radius, defined as the radius at which 
the average halo density is A c times the background den- 
sity according to the spherical top-hat model, is shown 
in column (5). A c is a number that depends on epoch 
and cosmological parameters (fi m ,fiA); for a flat ACDM 
model, A c <~ 337 at z = 0, where z is the redshift. The 
maximum circular velocity defined as 



^ GM(< r) 



1/2 



(5) 



where G is the gravitational constant and M(< r) is the 
mass contained within the radius r, is placed in column 
(6). In column (7) we show the C1/5 halo concentration 
parameter (see §4.1). The spin parameter A' as defined by 
Bullock ct al. (2001) is shown in column (8) (see §4.2). In 
column (9) we introduce the halo name tag for each one of 
the simulated halos. Cluster- and galaxy-sized halos are 
named Cl x and Gn x , respectively, where x denotes the 
statistics and n the halo number. 

4. RESULTS 

Figure 1 shows projected particle positions at z — from 
the 100h^ 1 Mpc box low-resolution simulations in slices of 
one-tenth of the box size. The differences in the parti- 
cle distribution between the lognormal realizations and 



the Gaussian ones are clearly seen: the density distribu- 
tion in the LN p statistics is lumpier with concentrated 
isolated peaks and a smaller coherence length, while in 
the LN n statistics it has a cellular structure with bigger 
voids and a larger coherence length. Notice that the void 
population is different in both cases; for the LN p model 
it appears to be larger than in the Gaussian case, while 
the contrary happens in the LN„ model. The contrasts 
among the different particle distributions are not as re- 
markable as shown, for example, in MMLM because of 
the different normalizations: we normalize our models to 
the same a% while MMLM stop their simulations when 
"the particle two-point function is best fitted by a power 
law £(r) = (r/r )~ 7 , with 7 = 1.8, in a suitable inter- 
val"; accordingly, their LN p (LN„) models are less (more) 
evolved. We measured £(r) at z — for our different statis- 
tics and found, as expected, similar correlation length ro 
(defined as the radius r where £(r) = 1) values. However, 
for r < r the correlation functions of the different sce- 
narios start deviating from each other, with the greatest 
differences found at the smallest scales. As expected, the 
£ functions of the LN p , G p (or G„) and LN n statistics are 
on the top, middle, and bottom, of the £ vs. r diagram, 
respectively. In any case, it is not the aim of this paper to 
explore the large-scale structure produced in simulations 
with NGICs, but to explore the effects of NGICs on the 
properties of CDM halo. 

Following, we present a study of the structural prop- 
erties of galaxy- and cluster-sized halos from the non- 
Gaussian realizations and their comparison with the cor- 
responding halos of the Gaussian simulations. 

4.1. Concentrations and Density Profiles of Halos 

In Figure 2 we show several general concentration 
parameters (vs. the halo mass M/j) measured for 
the halos in the low resolution simulations (12.5, 20 
and 100 h _1 Mpc box size). These plots give us a 
preliminar idea of the effects of NGICs on the halo 
properties for a relatively large number of objects. 
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Fig. 2. — Different concentration parameters (see text) vs. 
M» for the halos in the low-resolution simulations (small sym- 
bols) and those re-simulated with high-resolution (large sym- 
bols). Only halos with more than 1500 particles are plot- 
ted. The different symbols indicate the box size and the 
statistics of the simulations (see the legends in the upper 
panel). The dashed lines in the upper and medium pan- 
els are the linear fitting to the data (Gaussian isolated ha- 
los) presented in Avila-Reese et al. (1999), and the analyti- 
cal estimate of cm assuming the NFW profile and using log 
Cvi r = 2.35 — 0.1log(M v /h~ Mq). According to each one of 
the three concentration parameters used, the LN P (solid sym- 
bols) and LN„ (open symbols) halos are systematically more 
concentrated and less concentrated than the Gaussian halos 
(squeletal symbols), respectively. 



Only halos with more than 1500 particles are plotted. The 
Ci/5 concentration (upper panel) is defined as the ratio be- 
tween the halo radius (virial or truncation radius 1 ) and 
the radius where is contained one-fifth of the halo mass 
(Avila-Reese et al. 1999). This parameter depends on 
the definition and determination of the halo radius and 
mass. It is convenient to introduce other concentration 



parameters that do not depend on these definitions. Here 
we use the Cm and cs parameters. The cm concentration 
(medium panel) is defined as 27 times the ratio between 
the mass at r; n = 8.5 kpc (v max /220kms~ 1 ) and the mass 
at r out = 3xri„ (Dave et al. 2001). The cs concentration 
(lower panel) is defined as the ratio between the radius 
where the halo mean over density is 3 x 10 4 and where it is 
800. The faster the density between these two radii decays 
(larger concentration), the larger is cs, tending to 1 in the 
limit. Notice that cm traces the halo concentration in 
the inner region, while c$ refers to a large portion of the 
halo, it is more global than cm- We do not use the more 
common NFW concentration parameter, c v i r , because we 
do not know a priori whether the non-Gaussian halos will 
be described by the NFW density profile. 

It is clear from Figure 2 that halos produced in the 
LN p model are systematically more concentrated than the 
Gaussian ones, while the LN„ halos show an opposite 
trend. The LN p halos from the different box realizations 
are plotted with solid symbos, while those from the LN„ 
simulations are plotted with open symbols. The corre- 
sponding Gaussian halos are plotted with skeletal symbols. 
Dashed lines are the expected values for NFW CDM halos. 
In the case of C1/5 we take the results reported in Avila- 
Reese et al. (1999) for the isolated halos, while for cm we 
use its numerical relations with c v ; r , where for c v i r we took 
their average values in function of M„ as measured in the 
Gaussian ACDM simulation studied in Avila-Reese et al. 
(1999): logc vir = 2.35 - 0.10 log(M„/h- 1 M Q ) 

Now we turn to the rc-simulated halos with high- 
resolution. Their concentrations are also shown in Fig. 
2 with the corresponding large symbols. Figure 3 shows 
their density profiles. Upper and lower panels arc for 
the skew-positive and skew-negative models, respectively. 
Left and medium columns are for the high- and low- 
concentrated galaxy-sized halos, respectively, and right 
column is for the single cluster-sized halo. The LN-model 
halo profiles arc plotted with solid lines, while the cor- 
responding Gaussian-model halo profiles (calculated with 
the same seed) are plotted with dashed lines. The dot- 
ted lines are the NFW fitting to the Gaussian halos; the 
fittings are satisfactory, excepting at the very inner re- 
gions where the halo profiles are slightly steeper than the 
NFW profile, in particular for the halos in the lower pan- 
els. As said in §2.2, we selected two different galaxy-sized 
halos from the low-resolution simulation in order to cover 
the intrinsical (cosmological) dispersion on halo structures, 
which is mainly due to the dispersion on the halo mass 
aggregation histories, given on its own by the statistical 
nature of the primordial density fluctuation field. 

The results presented in Fig. 3 show that the NGICs 
have a clear imprint on the density profile of the ACDM 
halos. For the LN p model the halos are not only more 
concentrated than the corresponding Gaussian-model ha- 
los, but they are also significantly cuspier. The opposite 
applies for the LN n -model halos. In order to see in more 
detail the features of the inner density profiles of the halos, 
in Fig. 4 we show their slopes vs. the radius. Upper and 
lower panels are for the skew-positive and skew-negative 
statistics. In the left panels the two galaxy-sized halos 



1 If the halo density profile flattens or even increases before attaining K v , then the halo is truncated at the flattening radius. This happens 
only in some cases for subhalos or multiple halos. In most of the cases R„is attained. 
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Fig. 3. — Density profiles of halos re-simulated with high resolution. Upper and lower panels are for the skew-positive and 
skew-negative statistics. The first and second columns are for the high- and low-concentration galaxy-sized halos, respectively, 
while the third column is for the only simulated cluster-sized halo. Solid lines are for the non-Gaussian halos and dashed lines are 
for the corresponding Gaussian counterpart (simulated with the same random see). Dotted lines are the best NFW fitting to the 
Gaussian-halo density profiles. The LN P (LN„) halos are more (less) cuspy than the Gaussian counterpart halos. 



are plotted (thick and thin lines are for the high and low 
concentration halos, respectively), while the right panels 
are for the cluster halo. As in Fig. 3, solid lines are for 
the non-Gaussian halos and dashed lines are for the corre- 
sponding Gaussian halos. In all the cases, the slopes of the 
non-Gaussian halos systematically depart from the slopes 
of the corresponding Gaussian halos at radii smaller than 
3 — 10% the virial radius. The inner slopes of the LN p 
(LN„) halos are smaller (larger) than those of their corre- 
sponding Gaussian halos. Notice that the Gaussian halos 
corresponding to the LN p run resulted with shallower pro- 
files than the Gaussian halos corresponding to the LN„ 
run (this is part of the cosmic variance); the latter have 



inner slopes around -1.4. 

4.2. Spin Parameter and Angular Momentum 
Distribution 

In Figure 5 we show the probability distribution of the 
modified spin parameter, p(X'), for all halos from the 12.5 
h _1 Mpc box with the different statistics; the G„ and G p 
halos were joined into one denoted by G (see top-left 
panel). The parameter A' determines the global angular 
momentum of a halo and it is defined as in Bullock ct al. 
(2001): 

A' = , (6) 

V2M V V V R V 
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Fig. 4. — Density profile slopes in the inner parts of the non-Gaussian and Gaussian halos showed in Fig. 3. Upper and lower 
panels are for the skew-positive and skew-negative statistics. The two galaxy-sized halos are shown in the right panels (thick and 
thin lines for the high- and low-concentrated halo, respectively) , while the cluster-sized halo is in the right panels. The line coding 
is the same than in Fig. 3. 



where J v is the angular momentum inside virial radius R^ , 
and Vv is the circular velocity at radius K v . In Table 1, 
column (7), we show A' for the high- resolution halos. The 
halo angular momentum is calculated as: 



i=l 



n x Vi, 



(7) 



where Tj and Vi are the position and velocity of the ith 
particle with respect to the halo center of mass. 

The spin parameter distributions shown in Fig. 5 were 
obtained using halos with more than 100 particles. An 
extra simulation with 256 3 particles, with the same pa- 
rameters as the 12.5 h _1 Mpc box simulation described in 
§3, was also run and stopped at z = 1. The about 1000 
halos with more than 500 particles from this simulation 
were used to compute a well resolved LN p spin-parameter 



distribution (top-right panel) . A complete analysis of this 
simulation has been deferred to a future paper. For Gaus- 
sian initial conditions, iV— body simulations have shown 
previously that p(X') can be well described by the lognor- 
mal distribution (e.g., Bullock et al. 2001), 



P(A') 



1 



er A 'V27r 



exp 



in 2 

2<j x , 



(8) 



We confirm this for our Gaussian simulations and find 
that the halo spin parameter distribution in a lognor- 
mal statistics (with positive or negative skewness) is also 
roughly lognormal. Curves on panels are the lognormal 
best-fits to the data: {ay, A[,)=(0.63,0.035), (0.61,0.029), 
(0.73,0.026), and (0.49,0.042) for Gaussian, LN p (high res- 
olution), and LN p and LN„ (low resolution) simulations, 
respectively. The figure shows that halos formed in a LN p 
statistics are biased to lower values of A'; on the contrary, 
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Fig. 5. — Probability distribution (histogram) of the halo spin parameter A' for different statistics: Gaussian (top-left panel), 
skew-positive and negative lognormal from the low-resolution simulations (bottom panels), and skew-positive lognormal from a 
high- resolution simulation stopped at z = 1 (top-right panel). Curves on each panel are lognormal best-fits to the data. The Ao 
parameter (see eq. [8]) in each curve is: 0.035 for the Gaussian (Gp and Gn added), 0.026 for LN P (0.029 for the high resolution 
z = 1 simulation), and 0.042 for LN n (the ax parameter oscillates around 0.6, see text). As can be seen from the plot and confirmed 
by the Aq values obtained from the fits, we find that LN n (LN P ) halos have on average a higher (lower) A' value than their Gaussian 
counterparts. 



halos in the LN„ statistics seem to have on average a higher 
A' value than their Gaussian counterparts. 

For the halos re-simulated with high resolution, we 
compute also the halo angular momentum distribution 
(AMD). Figure 6 shows a comparison between the AMD 
of the galaxy-sized halos simulated with the skew-positive 
statistics, LN p (left panels) and skew-negative, LN„ (right 
panels), and the AMD of each one of their Gaussian coun- 
terparts. The AMD is computed as follows. We find first, 
for each particle, the specific angular momentum compo- 
nent along J Vl calling it j. We then divide the sphere of 
radius R v in spherical shells and each of these shells is in 



turn divided in four sections. Differences in the number of 
particles between cells for most of them are below a factor 
of two. The total j of any cell can be negative because 
j is a projected component. For each halo, we measured 
the mass fraction that there is in cells with negative j and 
found that it ranges from 5% to 50%. It is remarkable 
to find a halo (halo G2q p ) that, with the cell geometry 
adopted, happens to have about half of its cells pointing 
to the opposite direction of its total angular momentum. 
The calculation of the AMD proceeds only for those halos 
for which this fraction is small, less than ~10%. Cells are 
ranked according to their j value (cells with negative j 
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Fig. 6. — Specific angular momentum j distribution of halos Gl y G2 for the two lognormal statistics LN P and LN„, and their 
corresponding Gaussian realizations. Only cells with positive j were considered and the total mass Mt was redefined to be the mass 
contained only in those cells with positive j. The line is taken from the analytical fit of Bullock et al. (2001) m(< j) — fij/1 + j, 
where m(< j) = M(< j)/M vir and j = j/jmax, with /i = 1.25. No systematic differences are seen between halos simulated with 
lognormal NGICs and those simulated with the Gaussian statistics. 



are rejected) and M(< j) profiles are then built by count- 
ing the cumulative mass in cells with angular momentum 
smaller than j. Mass and angular momentum are given in 
M„ and j max units, respectively, where j ma x is the maxi- 
mum value reached by j. We have also plotted in Figure 
6 the analytical form of M(< j) proposed by Bullock et 
al. (2001) with ji = 1.25; the line is intended to guide the 
eye, it is not a fit to our numerical results. 

4.3. Ellipticities 

We have measured the ellipticities of halos using the 
tensor of inertia. This is defined as 

hi = x^j/r 2 , (9) 

where the sum is over all particles within r vir , Xi (i = 
1,2,3) are the particle coordinates with respect to the halo 



center of mass, and r is the distance of the particle to the 
halo center. The ellipticities are then given by 

ei = l-— , e 2 = l-— , (10) 

^3 ^3 

where A3 > A 2 > Ai are the eigenvalues of the tensor of 
inertia. We evaluate the triaxiality parameter using the 
following formula (e.g., Franx, Illingworth, & de Zeeuw 
1990) 

T=^ ^ (11) 

A halo is prolate (oblate) if T = 1.0 (T = 0.0). 

The ellipticities e\ (solid lines) and e 2 (dashed lines) as 
a function of radius are shown in Figure 7 for halos G2Q n 
and G2lnp (left panel) and for halos C1g p and CIlnp (right 
panel). For these halos, no systematic difference is seen in 
ellipticities when using different statistics. 
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Fig. 7. — Ellipticities el and e2 as a function of radius are shown in left (halo G2) and right (halo CI) panels. By definition 
el < e2 and thus el curves (dashed lines) lie below the corresponding e2 curves (solid lines). 



4.4. Subhalo Population 

In Figure 8 we show the present-day subhalo v max cu- 
mulative function for all our galaxy-sized halos displayed 
in Table 1. All subhalos inside a radius of 200 h _1 kpc 
from the center of their parent halo are counted. We plot 
this function down to v max ~ 20 kms -1 , where we believe 
we are complete (Klypin et al. 1999). The observational 
data (five pointed stars) are taken from Klypin et al. Two 
interesting and somewhat expected results come just from 
looking at the figure: substructure is suppressed in ha- 
los simulated with a LN p statistics, as compared with the 
Gaussian one, while a slightly higher number of subha- 
los survive in the LN„ simulations. These results can be 
explained as an effect of the stronger (lighter) tidal field 
developed by their parent halos in the corresponding LN p 
(LN„) statistics (see §3.1 and Colin, Avila-Reese, & Valen- 
zuela 2000). 

5. SUMMARIZING DISCUSSION 

We have analyzed the structure and substructure of 
ACDM halos obtained in low and high-resolution cos- 
mological A— body simulations using non-Gaussian ini- 
tial conditions (NGICs). The NGICs were generated in 
the gravitational potential with a lognormal statistics, 
positive-skewed (LN P ) in one case and negative-skewed 
(LN„) in the other. The sign of the skewness remains the 
same when passing to the density PDF but the statistics 
are no longer lognormal. Since our aim is to explore in a 
generic way the influence of the NGICs on the properties 
of the halos, the particular choice of the NGICs is not rele- 
vant. In order to attain a fair comparison with results from 
the usual Gaussian-initial-conditions simulations, for each 
non-Gaussian simulation, we simulated the corresponding 
Gaussian one, using the same random seed. After the 
Zel'dovich approximation (z = 50), the skewnesses of the 
Gaussian, LN p , and LN„ PDFs measured in the simula- 



tions are 0.52, 3.90, and -1.27, respectively, while the G3 
parameter values are 3.05 and 0.067 for the LN p and LN„ 
cases, respectively (for the Gaussian case G3 = 1 by defini- 
tion, see eq. [1]). These values are after some gravitational 
evolution, and they fall within the range of several of the 
theoretically proposed non-Gaussian statistics. 

As discussed in the Introduction, the question of Gaus- 
sianity or not in the primordial fluctuation field is highly 
debated in theory as well as from the point of view of 
the observations. Most of the evidence shows that if the 
primordial fluctuation field is not-Gaussian, then the devi- 
ations from Gaussianity should not be very large. On the 
other hand, there is no reason to a priori conclude that the 
statistics is the same at all the scales. Most of the tests 
on Gaussianity were applied for relatively large scales (in 
the CMBR anisotropy maps and for cluster abundances). 
Even in this range of large scales, there have been some 
observational hints that the statistics could be different 
at different scales (e.g., Chiang et al. 2003). As far as 
we know, the present paper introduces for the first time 
a potential way to test Gaussianity at the smallest scales, 
using for this a comparison of the structure of galaxy and 
cluster-simulated halos with observational inferences. 

Our main conclusion is that the structure and substruc- 
ture of the galaxy- and cluster-sized ACDM halos are 
affected by the initial statistics. For the skew-positive 
(skew-negative) statistics, halos are more (less) concen- 
trated, cuspier (shallower), and with less (more) substruc- 
ture than the corresponding Gaussian halo counterparts. 
These main findings are clearly seen in Figs. 2 and 3. 

According to our results, the skew-positive NGICs 
would exacerbate the difficulty that the Gaussian ACDM 
model faces with observations: the rotation curves of 
dark-matter dominated (dwarf and low surface bright- 
ness) galaxies show that the inner density profiles of the 
halos are shallower and less concentrated than those of 
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Fig. 8. — Cumulative subhalo v max — functions of halos Gl y G2 for the Gaussian and lognormal statistics. All subhalos with 
v max greater than 20 kms -1 within a sphere of radius 200 h _1 kpc were counted. The averaged v max — function from satellites of 
the Milky Way and Andromeda is represented by stars (taken from Klypin et al. 1999). We plot for clarity only the error bars for 
the "observed" circular velocity function. It is clear that substructure is erased more efficiently in halos with a LN P statistics than 
with a Gaussian one (see top panels). This effect can be traced back to the stronger tidal field experimented by subhalos in the 
LN P scenario. The opposite effect is expected when the amount of substructure in a LN„ scenario is compared with its Gaussian 
counterpart. This is indeed the case: more subhalos survive under the LN n statistics than under the Gaussian one. 



the ACDM halos (see Colin et al. 2002 for references). 
We notice that tests related to cluster abundances and 
their evolution suggested skew-positive deviations from 
Gaussianity (see Introduction). On the other hand, the 
skew-negative NGICs work in the direction of making the 



galaxy-sized dark halos less concentrated and shallower 
than the Gaussian ones. 

Regarding the subhalo population, since the halos in 
the LN p statistics are very concentrated, more subhalos 
are tidally destroyed than in the Gaussian case. The op- 
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posite occurs with the subhalo population in the LN„ sim- 
ulations, although the difference with the Gaussian re- 
sults in this case is very small. Observations show that 
the Gaussian ACDM model overpredicts the number of 
satellite subhalos within simulated galaxy-sized CDM ha- 
los (Klypin et al. 1999; Moore et al. 1999a). Neverthe- 
less, this difference could be due to astrophysical processes 
(reionization, for example) which inhibit the formation of 
dwarf galaxies within the subhalos (Bullock, Kravtsov, & 
Weinberg 2000, Benson et al. 2002). 

The halo spin parameter probability distribution also 
has shown to depend on the primordial statistics: LN„ 
(LN p ) halos have on average a higher (lower) A' value than 
their Gaussian counterparts. A potential difficulty pointed 
out for the Gaussian halos is that the disks formed within 
them result too much concentrated (Navarro & White 
1994). A possible way to improve on this problem is to 
generate halos with larger spin parameters than the Gaus- 
sian CDM ones. For example, this is the case for the LN„ 
model presented here. 

From a theoretical point of view, our results raise an 
interesting problem: the density profile of the dark matter 
halos seems to depend on initial conditions. The question 
whether the structure of CDM halos is universal and non- 
dependent on cosmology, power spectrum, formation his- 
tory, etc. has been several times discussed in the literature. 
Avila-Reese et al. (1998, 1999) have shown that the CDM 
halo density profiles oscillate around the NFW profile, de- 
pending mainly on the halo mass aggregation history (see 
also Wechsler et al. 2002). Nevertheless, a significant 
variation of the inner density profile was not expected, 
unless thermal processes are included, as is the case of 
self- interacting dark matter (e.g., Colin et al. 2002) or 
high rms- velocity dispersions due to thermal relicts (e.g., 
Avila-Reese et al. 1998, 2001; Lokas & Hoffmann 2000). 
It has been also shown that variations in the power spec- 
trum shape and the type of the halo collapse (hierarchical 
merging or by fragmentation) do not introduce significant 
changes in the halo density profiles (Moore et al 1999b; 
Huss, Jain, Steinmetz 1999; Avila-Reese et al. 2001; Man- 
riquc et al. 2003). The fact that the initial statistics 



of the fluctuation field does significantly affect the inner 
structure of the CDM halos is somewhat unexpected al- 
though, using the semi-numerical method of Avila-Reese 
et al. (1998), it was already shown that for a possitive log- 
normal statistics (in the density field), the halos are more 
cuspy than in the Gaussian case (Viniegra 1997; Avila- 
Reese 1998). Halos form typically from the high-density 
peaks in the matter distribution. In the skew-possitive 
statisitics, high-density peaks are more abundant, isolated 
and "coherent" than in the Gaussian statistics; the oppo- 
site happens with the skew-negative statistics. The final 
structure of the halos is probably related to these facts. 
A deeper analysis of this question has been deferred to a 
future paper. 

In the last few years, the observational determination of 
the cosmological parameters attained an unprecedent high 
accuracy and the body of astronomical observations of the 
local and high-redshift universe has grown dramatically. 
This progress not only allows to constrain better the cos- 
mological models, but also opens a new dimension in the 
relaxation of some hypotheses made before, as well as in 
the exploration of second-order phenomena in the cosmic 
structure formation process. An example of this is just the 
statistics of the initial fluctuation field. Having now well 
defined the main cosmological parameters and a succcsful 
model of cosmic structure formation, one may explore the 
influence of NGICs on structure formation free of degen- 
eracies with other parameters and assumptions (e.g., with 
Q m and cr 8 ). 
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